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Abstract. Theoretical spectral energy distributions are 
computed for a grid of hydrogen-deficient and carbon-rich 
model atmospheres with a range in T e g of 4000 — 5500 K 
and log g of 1.0 — 0.0, using the technique of opacity sam- 
pling, and taking into account continuous, atomic line 
and molecular band absorption. The energy distributions 
are compared with the spectrum of V4334 Sgr (Saku- 
rai's object) of August 12, 1998 in the wavelength inter- 
val 300 — 1000 nm. This comparison yields an effective 
temperature of V4334 Sgr of T cff = 5250 ± 200 K, and a 
value for the interstellar plus circumstellar reddening of 
£ B -v = 1.3 ±0.1. 
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1. Introduction 

V4334 Sgr, the "novalike object in Sagittarius", was dis- 
covered by Y. Sakurai on February 20, 1996 (Nakano et 
al. 1996). Soon after discovery, it was found to be a final 
He flash object, a star on the rare evolutionary track that 
leads from the central stars of planetary nebulae back to 
the red giant region, an object also called a "born-again 
giant" . Its progenitor was a faint blue star (~ 21 m ) in the 
centre of a low surface brightness planetary nebula (Duer- 
beck & Benetti 1996). Early spectroscopic studies of the 
object found an increasing hydrogen deficit in the atmo- 
sphere, and a C/O ratio > 1 (Asplund et al. 1997, Kipper 
& Klochkova 1997). 

Unlike the slowly evolving born-again giant FG Sge, 
which rose in brightness during the first half of the 20th 
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Fig. 1. The visual light curve of V4334 Sgr, according to 
Duerbeck et al. (2000). The times of five spectroscopic 
observations, analyzed by Asplund et al. (1997, 1999) and 
Kipper & Klochkova (1997), are indicated by open trian- 
gles. The times of our two spectroscopic observations of 
April, 1997 and August, 1998 are marked by filled trian- 
gles. The second observation took place at the beginning 
of a deep brightness decline. 



century, V4334 Sgr is a very quickly evolving object 
(Fig. 1). It shows striking similarities with the unfortu- 
nately poorly documented final He flash object V605 Aql 
that erupted in 1918 - 1923 (see Clayton & de Marco 
(1997) for an overview). Between March 1996 and April 
1997, V4334 Sgr evolved from early F to late red giant 
stages. Infrared observations by Kimeswenger et al. (1997) 
and by Kamath & Ashok (1999) show a noticeable excess 
in the infrared K band starting March 1997, indicative 
of dust formation. However, only in early 1998 a ~ l m 
depression appeared in the optical light curve (Liller et 
al. 1998, Duerbeck et al. 1999, 2000). This moderate de- 
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pression persisted for a few months, and was followed by 
a ~ 6 m decline which started in August 1998. After a 
partial recovery in late 1998 - early 1999, the dust obscu- 
ration continued with increasing strength, and it may be 
assumed that in the year 2000, only the short-wavelength 
tail of the cool dust envelope can be observed in the opti- 
cal spectrum, while the stellar photosphere is engulfed in 
the opaque dust shell. Unfortunately, not much is know 
about these late phases, since the spectroscopic coverage 
is poor due to the extreme faintness of the object in the 
visible region. 

A spectrum of V4334 Sgr was obtained on April 29, 
1997, before any dust obscuration in the optical region be- 
came evident, but after the strong IR excess had evolved. 
This spectrum was compared with theoretical spectral en- 
ergy distributions (Pavlenko, Yakovina & Duerbeck 2000). 
The authors derived a temperature of T e g = 5500 K, and 
an interstellar reddening of £b-v = 0.7. In this paper, we 
apply the same methods to a spectrum which was taken 
15 months later. 

2. Observations 

On August 12, 1998, a spectrum (range 355 — 1000 nm, 
resolving power ~ 2000) was taken by A. Piemonte with 
the B&C Cassegrain spectrograph attached to the ESO 
1.52 m telescope, and kindly put at the authors' disposal. 
It was taken around the time of the onset of the ~ 6 m 
decline, at a time when V4334 Sgr was obviously already 
suffering ~ l m of visual dust obscuration (Fig. 1). 

The spectrum was reduced using standard IRAF rou- 
tines, and the spectral energy distribution (SED) was de- 
rived using spectrophotometric standard stars. The ob- 
served SED is influenced by interstellar reddening as well 
as by additional circumstellar reddening. The value of the 
interstellar reddening is still not accurately known; we will 
assume £b-v = 0.7 in the following, as was derived in the 
analysis of the "dustless" spectrum of April 29, 1997, by 
Pavlenko, Yakovina & Duerbeck (2000); note, however, 
that Duerbeck et al. (2000) prefer £b-v — 0.8, a value 
which is also consistent with the observed versus theoreti- 
cal Balmer decrement of the planetary nebula surrounding 
V4334 Sgr (Kerber et al. 2000). 

We will assume that the dusty envelope acts like an 
additional amount of interstellar reddening. Dcreddened 
spectral energy distributions of the August 12, 1998 spec- 
trum were calculated for assumed values 0.7 < £b-v < 
1.3, and compared with model atmospheres. 

3. Procedure and Results 

Our computations of model atmospheres and synthetic 
spectra are carried out in a classical approach: a plane- 
parallel model atmosphere in LTE, with no energy diver- 
gence. We assumed chemical compositions of V4334 Sgr 
as determined by Asplund et al. (1997) and Kipper & 



Table 1. Abundances of H, He, C, N, O (scaled to J2 N i 
= 1) used in this paper. 



Element 


Asplund 


Asplund 


Kipper & 


[C] = +0.6 




et al. 


et al. 


Klochkova 






(1997) 


(1999) 


(1997) 




H 


-1.730 


-2.42 


-2.10 


-2.45 


He 


-0.027 


-0.020 


-0.01 


-0.05 


C 


-1.73 


-1.62 


-2.01 


-1.05 


N 


-2.53 


-2.52 


-2.70 


-2.55 


O 


-1.93 


-2.02 


-2.59 


-2.05 



Klochkova (1997). Their results differ in some details, last 
not least because there are differences in the physical in- 
put parameters and in the details of the procedure. 

The abundances of some elements obtained by As- 
plund et al. (1997) for May 1996 are given in the sec- 
ond column of Table 1. Asplund et al. (1999) obtained 
another set of abundances for October 1996. For com- 
parison, we show the abundances of Kipper & Klochkova 
(1997), obtained independently for a spectrum of July 
1996. In general, their results differ more from Asplund 
et al. (1997) than the results of Asplund et al. (1997) and 
(1999). Therefore, they provide a more interesting base- 
line for numerical experiments carried out to clarify the 
dependence of our results on input abundances. 

From the common point of view, hydrogen is an im- 
portant clement, and its abundance in the atmosphere 
of V4334 Sgr is steadily declining (Asplund at al. 1997). 
However, Pavlenko et al. (2000) showed that the depen- 
dence of our results on log iV(H) is much weaker than that 
found by Asplund et al. (1997) for the early stages of its 
evoluton. Indeed, in the case of lower T e g (< 6500 K) the 
contribution of H and H _ to the total opacity in the line 
forming region is reduced in comparison with "hot mod- 
els" , because molecular absorption which does not include 
hydrogen compounds dominates there. 

Furthermore, and more importantly, the input carbon 
abundance of Asplund et al. (1997) is larger by up to 0.6 
dex than that which was obtained by them from the fine 
analysis of high resolution spectra. That is the well known 
"carbon abundance problem" (Asplund et al. 2000), which 
has, unfortunately, not yet been resolved. 

One may think that, in the meantime, the model at- 
mospheres should not be computed with the derived C 
abundances by Asplund et al. (1999) but rather with the 
higher input values. The use of the derived lower C abun- 
dance would result in a new "carbon problem" and thus 
even lower C abundances would be derived. To clarify the 
case, we studied the impact of higher carbon abundances 
on the model atmospheres and SEDs of V4334 Sgr. We 
used 4 abundance sets (i.e. Asplund et al. (1997) for May 
1996, Kipper and Klochkova (1997) for July 1996, Asplund 
et al. (1999) for October 1996, and the "carbon rich case" 
with [C] = + 0.6 with respect to the Asplund et al. (1999) 
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abundance) in order to study different aspects of the prob- 
lem: 

— to model the impact of abundance changes on our re- 
sults, 

— to avoid possible uncertainties related with the "car- 
bon problem" . 

Opacity sampling model atmospheres and syntheti- 
cal spectra of V4334 Sgr were computed by the SAM941 
(Pavlenko 1999) and WITA6 programmes, respectively. 
The model atmosphere computation procedure is de- 
scribed in Pavlenko et al. (2000) and Pavlenko (2000). To 
clarify the topic, some details of the computations proce- 
dure are given here: 

— Convection was treated in the frame of the mixing- 
length theory with l/H = 1.6. Overshooting was not 
considered in the computations; 

— the list of atomic lines was taken from VALD 
(Piskunov et al. 1995); 

— molecular opacities were computed in the frame of 
JOLA approach. We take into account the absorption 
in the 20 band system of diatomic molecules (see Ta- 
ble 2 in Pavlenko et al. 2000); 

— bound-free absorption of C I, N I and O I atoms were 
computed using cross-sections of Hoffsacss (1979); 

— SAM941 and WITA6 use the same opacity lists. 

Model atmospheres, synthetical spectra and SEDs 
were computed for a microturbulent velocity of 5 km/s, 
which is a typical value for atmospheres of post-AGB 
stars. The resulting theoretical spectra were convolved by 
a gaussian with a half- width of 0.5 nm. 

Basically, the atmospheric structure of chemically 
peculiar stars should respond to abundance changes 
(Pavlenko & Yakovina 1994, 2000). To study the impact 
of the set of input abundances on our results, we carried 
out a few numerical experiments with model atmospheres 
computed with abundances of Asplund et al. (1999) and 
Kipper & Klochkova (1997). 

The comparison of the temperature structure of a 
model atmosphere T e s/\ogg — 5200/0.0 computed for 
the four abundance grids is shown in Fig. 2. We prefer to 
use pressure as a depth parameter, because temperature 
structures in the coordinates (r ross , T) look very similar 
(see Fig. 2b in Pavlenko et al. 2000). In general, the photo- 
spheres are shifted outwards for higher carbon abundances 
due to the opacity increase. Let us note a few items: 

— The structure of model atmospheres of V4334 Sgr, 
computed for abundances by Asplund et al. (1999), 
Kipper & Klochkova (1997), and the "carbon- rich 
case" differ substantially, especially in the deep lay- 
ers. Differences of model atmospheres computed for 
Asplund et al.'s (1999) abundances obtained for May 
and October of 1996 are rather marginal; 
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Fig. 2. The T = f(P g ) structure of model atmospheres 
computed for Asplund et al. (1997), Kipper & Klochkova 
(1997), Asplund ct al. (1999), and [C] = + 0.6 abundances. 
A second diagram shows their T = /(r) structure. 

— the structure of the outermost layers of model atmo- 
spheres computed for the Asplund ct al. (1997) and 
Kipper & Klochkova abundances is similar. Moreover, 
the differences of the computed SEDs are rather weak 
(Fig. 3). More pronounced differences are found for the 
"carbon-rich" case; 

— in the blue region, the "flux peaks" which correspond 
to the minima of molecular absorption, are less pro- 
nounced for the model based on Kipper & Klochkova's 
(1997) abundances, but even more so in the observed 
spectrum (Fig. 3); 

— We have to admit that there is an obvious lack of com- 
puted opacity in the blue region, and the model is only 
correct in a qualitative sense (Pavlenko & Yakovina 
2000). 
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Fig. 3. SEDs computed for model atmospheres 5200/0.0 
with abundances taken from Asplund et al. (1999), Kipper 
& Klochkova (1997) and the "carbon-rich" case. The first 
two SEDs nearly coincide. 
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Fig. 4. Synthetical spectra for different effective temper- 
atures and logg = are superimposed on the spectrum of 
V4334 Sgr (thick line) observed in August, 1998. 



The overall shape of the spectrum of V4334 Sgr is gov- 
erned by the bands of the C2 Swan system and by the 
violet and red systems of CN (Pavlenko et al. 2000). In 
the near UV (A < 400 nm), atomic absorption becomes 
important (see Pavlenko & Yakovina 2000). 

Computed and observed spectra were normalized to 
equal flux at A ~ 740 nm. In general, our main conclu- 
sions show a rather weak dependence on the choice of A 
normalization (see Pavlenko et al. 2000). 



Nevertheless, the SED of V4334 Sgr depends critically 
on T ff (Pavlenko & Yakovina 2000). For the spectrum 
of August 12, 1998, a best fit was obtained for T e g ~ 
5250 K and E B -v = 1.3 (Fig. 4). The fit is based on the 
Asplund et al. (1999) abundances for October 1996. The 
corresponding values derived for our spectrum of April 
29, 1997, are T cff = 5500 K and £ B -v = 0.7 (Pavlenko 
et al. 2000). In contrast to the rapid cooling of V4334 Sgr 
in 1996, T e g has apparently declined by only ~ 250 K 
in the course of the 15 1 months between early 1997 and 
mid-1998. 



4. Discussion 

Our SEDs were computed in the plane-parallel approach. 
This may be used for T c g — 5500 K (sec discussions in 
Asplund et al. 1998 and Pavlenko et al. 2000), but model 
atmospheres for T c s = 5250 K are more likely affected by 
sphericity effects. Fits of model atmospheres with logg = 
and log g = 1 to the observed spectrum have a similar 
quality. This, however, only indicates that the dependence 
of the model on gravity is rather weak. 

It is noteworthy that our theoretical fluxes are much 
too high between 600 and 700 nm, which suggests that 
an important molecular line opacity is missing. A similar 
effect was noted by Asplund et al. (1997) in R CrB for A 
500 nm. It is doubtful that the assumption of geometry 
and homogeneity is able to produce a shortage of flux at 
such a very specific wavelength interval. In V4334 Sgr, the 
effect is probably caused by the absorption of one or more 
molecules which are not yet identified. 

The observed and computed spectra for August, 1998 
are more discrepant in the blue region, while the fit in the 
red looks acceptable. A similar phenomenon was found in 
the investigation of the April 1997 spectrum. By August 
1998 the dust envelope had become thick, while it was ab- 
sent or marginal in April 1997. Thus we conclude that the 
discrepancies are likely not caused by dust, but by one 
of the other effects outlined above. In the frame of our 
simple model, which assumes that the circumstellar and 
the interstellar dust have similar properties in the visible 
part of the spectrum, £b-v increased from 0.7 to 1.3, ac- 
cording to our spectroscopic analysis. Apparently, in 1997 
there was no "active" optical circumstellar reddening in 
the line of sight, which is proven by the similar effective 
temperatures derived from spectroscopy and photometry 
(Fig. 2 of Asplund et al. 1997), while at that time some 
dust clouds may have already formed outside the line of 
sight, in order to account for the IR excess as observed by 
Kimeswenger et al. (1997) and Kamath & Ashok (1999). 
The increase in £b-v between April 1997 and August 
1998, A_Eb-v = 0.6, is thus caused by circumstellar red- 
dening of newly formed dust. If the optical properties of 
this material are similar to those of interstellar material, 
it should cause a drop in V by 0.6 x 3.1 ~ 1.85 mag 
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Fig. 5. Fits to SEDs with model 4000/0.0 of different 
abundances. The expected very strong carbon bands can 
be decreased in strength by reducing the C/O ratio, but 
the overall fit would be still poor. 

(with an uncertainty of ±0.1) with respect to a "dust- 
free" lightcurve. 

In order to estimate the magnitude of V4334 Sgr in 
such a "dustfree" lightcurve, the trend in the 1996-97 V 
light curve (as taken from Duerbeck et al. 2000) was lin- 
early extrapolated to August, 1998, and an "expected" 
magnitude of V = 10.75 was derived; if only the 1997 data 
are used for an extrapolation, the result is V — 10.65. In 
fact, the object is observed at V = 12.55. If we trust in 
the validity of the linear extrapolation, the circumstcllar 
Ay, which is necessary to put the object to the observed 
brightness, is 1.85±0.05, in good agreement with the value 
expected from the result of the spectroscopic study. 

Is there another way to explain the observed spectrum 
of August, 1998, e.g. by keeping the reddening at its orig- 
inal value of £b-v — 0-7, and decreasing the surface tem- 
perature? If T c g is assumed to be as low as 4000 K, a 
general agreement with the overall SED may be achieved, 
but the fit to the molecular bands is extremely poor. The 
fit cannot be improved, even if the abundances of C, N, O 
are permitted to vary (Fig. 5). Here we show computations 
for different C/O ratios (given in logarithmic scales). If O 
is more abundant than C, a completely different molecu- 
lar chemistry, and therefore bands of oxygen compounds 
would appear in V4334 Sgr. Fig. 5 illustrates that any 
changes of carbon abundances do not result in a better fit 
of the observed spectrum. 

The estimate of T e g depends on the assumed -Eb-v- 
Fig. 6 shows some fits of computed spectra to the observed 
spectrum of V4334 Sgr, dereddened for E B -v = 1.0. There 
we also plot the SED of V4334 Sgr, computed for the 
"carbon-rich" case. This permits us to better assess the 
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Fig. 6. Fits to SEDs of August 12, 1998 dereddened for 
E B -v = 1.0 and 1.3 with model 5250/0.0 of the "carbon- 
rich" case. There are some obvious problems of the fittings 
in the red and in the blue. 
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Fig. 7. Temporal changes of T c g in V4334 Sgr. Open 
circles: Asplund ct al. (1997), open squares: Asplund et 
al. (1999), open diamond: Kipper and Klochkova (1997), 
filled circle: Pavlenko et al. (2000), filled square: this pa- 
per. 



uncertainties in these estimates. As we see in Fig. 6, our 
fit with £b-v — 1.0 looks less perfect both in the red and 
in the blue both for the Asplund at al. (1999) abundances 
for May 1996 and even for the "carbon-rich" case. Neither 
variations of logAT(C) nor of £b-v permit to obtain a 
better fit to the SED observed in August, 1998. 
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Thus we conclude that T c g of V4334 Sgr was near 
5250 K in August, 1998, and that the overall modifica- 
tion of the SED is caused by circumstellar dust. 

The temporal changes of T e g during 1996-1998 are 
shown in Fig. 7. Following the assumption of Duerbeck et 
al. (1998), that changes of T e g were caused by the growth 
of the pseudophotosphere surrounding a mass-losing ob- 
ject that radiates at approximately constant luminosity, 
we conclude that its expansion had noticeably slowed 
down or had even almost come to a halt during 1997 
- 1998, and its temperature had almost stabilized. This 
stage was accompanied by dust formation processes in the 
envelope above the stellar photosphere, whose strength in- 
creased with time. The object recovered only partly from 
these events, so that by the second half of 1999 the photo- 
sphere of V4334 Sgr had become basically invisible in the 
optical region. 

5. Conclusions 

A comparison of the observed spectra of V4334 Sgr in 
April, 1997 and August, 1998 with computed SEDs shows 
that the temperature of the radiating atmosphere had only 
decreased in T e g by 250 K, while the £b-v had increased 
by 0.6. This appears to be in perfect agreement with the 
observed depression in the visual light curve, if a standard 
interstellar extinction law Ay = 3.1 x £b-v is assumed for 
newly formed circumstellar dust. Thus there is not only 
photometric, but also spectroscopic evidence that (a) the 
light curve depressions are caused by dust, and that (b) 
the luminosity and mass loss changed in such a way that 
the pseudophotosphere was kept at an almost constant 
effective temperature in 1997 - 1998. 
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